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Abstract 

We find static solitons stabilized by quantum corrections in a (l+l)-dimen- 
sional model with a scalar field chirally coupled to fermions. This model does 
not support classical solitons. We compute the renormalized energy functional 
including one-loop quantum corrections. We carry out a variational search 
for a configuration that minimizes the energy functional. We find a nontrivial 
configuration with fermion number whose energy is lower than the same num- 
ber of free fermions quantized about the translationally invariant vacuum. In 
order to compute the quantum corrections for a given background field we 
use a phase-shift parameterization of the Casimir energy. We identify orders 
of the Born series for the phase shift with perturbative Feynman diagrams 
in order to renormalize the Casimir energy using perturbatively determined 
counterterms. Generalizing dimensional regularization, we demonstrate that 
this procedure yields a finite and unambiguous energy functional. 



1. INTRODUCTION 

Many quantum field theories contain solitons, nontrivial configurations that are local 
minima of the energy functional. Often, the existence of a soliton can be anticipated from 
topological properties of the field theory. In such cases, topologically nontrivial configura- 
tions cannot be continuously deformed into a vacuum configuration. Still, a detailed look at 
the dynamics is required to see that such configurations do not collapse to a point. For ex- 
ample, the soliton in the Skyrme model is unstable without the addition of a four- derivative 
term to the action. In general, the question of soliton stability requires a detailed study of 
the energy. If we have a reliable method for computing the leading quantum contributions 
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to the energy of nontrivial field configurations, we can begin to study the stability of topo- 
logical and nontopological solitons. The techniques introduced in Ref. |j]-f| provide such a 
framework in renormalizable quantum field theories. 

In this paper we apply these techniques to reexamine suggestions made in previous 
works |5|,[|] that a heavy fermion can create a soliton. We consider a simple renormalizable 
quantum field theory in 1+1 dimensions. In this model, a two-component scalar field couples 
to fermions in a chirally invariant way. The fermion mass arises from a Yukawa coupling to 
the scalar, which has a nonzero vacuum expectation value. At the classical level the theory 
has no solitons. (This is in contrast to iV = 2 supersymmetric models considered in Ref. [0 , 
which support classical solitons that are unmodified by quantum corrections.) We investigate 
whether configurations of the scalar field that carry fermion quantum numbers can be lighter 
than ordinary fermion states built on top of the translationally invariant vacuum. We show 
that they can, and that the lightest of these configurations can be identified as a stable 
soliton. 

It is well known that a fermion can be strongly bound in the background of a spatially 
varying scalar field. But the binding energy comes with the cost of the classical energy of 
the scalar configuration. In addition, there is the possible cost of an increase of the fermion 
Casimir energy, since the fermion vacuum is polarized by the background scalar field. We 
must include this contribution because it is of the same order in the loop-counting parameter 
as the binding energy of the strongly bound level. 

We give a detailed description of our method for calculating the one-loop fermion con- 
tribution to the energy of a fixed scalar background. Our method is exact for an arbitrary 
spatially varying background. It sums all orders in the derivative expansion. Since we work 
in a renormalizable theory, our results are unambiguous; we use the same counterterms in 
the presence of the soliton that we determined from the renormalization conditions in the 
perturbative sector. In order to combine the infinite loop energy and the infinite counterterm 
contribution into a finite answer, we must regularize the theory. Phase shifts are a key in- 
gredient in our approach. The Casimir energy is the sum over the energy shift of each mode 
in the fixed background. We express the sum over modes as an integral over the continuum 
density of states, which we relate in turn to the phase shifts. Then the ultraviolet- divergent 
pieces of the loop energy are linked to low orders of the Born series for the phase shifts, which 
can be unambiguously identified with Feynman diagrams. In this framework the quantum 
corrections to the energy are expressed as the sum of two finite pieces: a momentum integral 
involving Born-subtracted phase shifts, and a small number of perturbatively renormalized 
Feynman diagrams, which correspond to the Born subtractions ||T]-§|. 

The (l+l)-dimensional model we consider has many features in common with the linear 
cx-model in 3+1 dimensions: renormalizability, chirally symmetric scalar-fermion interac- 
tions, the possibility of a Dirac Hamiltonian with a spectrum that is asymmetric in energy, 
and a topological structure. Hence our findings in the (l+l)-dimensional case suggest that 
similar phenomena might occur in 3 + 1 dimensions: heavy fermions can be realized as 
soliton configurations carrying fermion number. This result is in agreement with decoupling 
theorems ||, which show that solitons supply the necessary fermion numbers to maintain 
anomaly cancellation after a fermion is decoupled by increasing its Yukawa coupling. 

In a theory with Np fermions, quantum corrections due to boson loops are suppressed by 
a factor of Np relative to the fermion quantum corrections. Therefore we work with Np large. 
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However, boson loops require special attention in a (l+l)-dimensional model with (classical) 
spontaneous symmetry breaking. Because of infrared singularities, the symmetry is actually 
restored by these corrections and it would be inappropriate to construct a soliton about a 
classical minimum. To circumvent this difficulty we introduce an explicit symmetry breaking 
term that is large enough to guarantee a unique vacuum state. Other treatments of solitons 
at the quantum level based on the inverse scattering method [p|-|T0|| cannot accommodate 
this essential ingredient. 

Our paper is organized as follows: In Section ^| we introduce the scalar sector of the model 
and explain the necessity of explicitly breaking the chiral symmetry. Section || contains the 
discussion of the fermionic sector. In Section [| we combine the scalar and fermion sectors to 
form the full expression for the energy. The corresponding numerical results are presented 
in Section [5[ We summarize our findings in Section ^. In Appendix [A] we justify the phase 
shift method using dimensional regularization. In Appendix [B] we demonstrate the behavior 
of the strongly bound fermion level in the WKB approximation. In Appendix |C| we explain 
why the soliton background need not be reflectionless, as claimed in previous works. Some 



of the results presented in this paper were summarized in Ref. [11 



2. THE BOSON SECTOR 

We consider a two component scalar field governed by the Lagrangian 

c B = \dJ-d^-v{$) (2.i) 

where <fi = (0i, 02)- The potential V contains the usual Higgs term and an explicit symmetry 
breaking piece proportional to 0i, 

V {<!>) = g(^" v + — ) "2 Of) ~ av (<l>i-v). (2.2) 

If a = 0, the U(l) transformation 

0i + # 2 — > e ilp (0x + i<p 2 ) (2.3) 

is a symmetry of Cb, but we take a big enough to avoid infrared singularities that would 
arise when this symmetry is spontaneously broken. For a > 0, V(4>) has a unique classical 
minimum at cl = (v,0). The modes fluctuating about the classical vacuum have masses 

m 2 = (A + a) v 2 and m 2 = av 2 . (2.4) 

Note that on the chiral circle, <p ■ <p = v 2 , this model reduces to the sine-Gordon model. 

Although we are principally interested in the stabilizing effects of fermions, we must first 
consider the infrared singularities in the boson sector that could invalidate our approach. 
We wish to study scalar field configurations that are deformations of the classical vacuum. 
However, it is well known that when a = 0, spontaneous symmetry breaking does not 
occur The absence of long-range order in 1+1 dimensions can be traced to infrared 

singularities associated with the propagator of the massless mode. Making a big enough 
gives a large enough mass to the would-be Goldstone mode to stabilize the classical vacuum. 
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The infrared problems can be seen in the renormalized one-loop effective potential, V e 



eff; 
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where the eigenvalues of the 2x2 matrix d 2 V / d<pid(f)j: 
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±0 2 -v 2 + 2av 2 /X) . 



(2.5) 



(2.6) 



We have renormalized V3 so that the vacuum expectation value at 



and the particle 



masses m a and m n remain unchanged. As a — > 0, m 2 — > 0, so V^f develops a logarithmic 
infrared divergence reflecting the fluctuations that restore U(l) invariance. But for any 
nonzero a, Vjg is well defined and has a unique minimum at <p c \. 

When we consider the fermion one-loop contribution, we find that when a becomes 
small the character of configurations that minimize the energy changes. To see this, note 
that V^g(0) = for any scalar field configuration on the chiral circle. As a — > 0, is forced 
to the chiral circle because any other choice has infinite energy. On the other hand, for a 
slowly varying that circles zero once as x goes from — oo to +oo, the classical energy goes 
to zero if stays on the chiral circle and the scale over which varies gets large. This 
nontrivial field configuration has zero boson contribution to its energy through one loop in 
the limit a — > 0. (Since the configuration varies slowly, integrating eq. Q2.5|) over x gives 
a good approximation to the one-loop effective energy; because we are in 1+1 dimensions, 
all higher derivative terms go to zero as the scale over which varies gets large.) When 
we consider the fermion one-loop contribution to the energy of this configuration, we find 
that it too contributes zero. Furthermore this configuration has fermion number one. Thus 
as a —>■ the model has a zero-energy, infinitely large soliton carrying fermion number. 
We believe this object is destroyed by the fluctuations that restore U(l) invariance. In our 
numerical search for solitons, we find a sharp transition as we vary a. At small a, we find 
that configurations with large widths are energetically favored, and the width is controlled 
by a. At large a we find configurations with moderate widths are favored, and the width is 
no longer sensitive to changes in a. We are interested in the latter. 



3. THE FERMION SECTOR 

The purpose of this section is fourfold. First, we set up the formal expression for the 
contribution of the fermion sector to the energy, which is expressed through the eigenvalues 
of the single-particle Dirac equation. Second, we set up the formalism for solving the Dirac 
equation in a general background <p(x). We show how to obtain the phase shifts, which 
characterize the continuum states. Third, we explain how the properly regulated and renor- 
malized fermion one-loop contribution to the energy can be computed in terms of the phase 
shifts and the bound state energies. And fourth, we explain how the fermion number of the 
background can be computed in terms of the phase shifts and bound states. 
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A. The Fermion Casimir Energy 

We introduce Np Dirac fermions coupled in a chirally invariant way to 0. We suppress 
the flavor label, but keep track of factors of Np as necessary, 

C F = % - [tt, *] 0i + i [f , 75^>] 2 ) (3.1) 

with \I/ = ^'7°. The fermions acquire mass m = Gv when the boson field takes the 
vacuum value c i. We have been careful about operator ordering in eq. ( |3. 1| ) because of 
its role in determining the fermionic contribution to the vacuum energy. The ordering of 
anticommuting fermion fields is fixed by requiring charge conjugation invariance. We define 
C0i = 0i, C02 = —02, and choose a Majorana basis for the Dirac matrices, 7 = <t 2 , 7 1 = 203, 
and 75 = Oi, so that C^> = \&*. With these definitions it is easy to verify that is even 

under C and that [\j/,7 5 \I/] and [\I/,7 M \I/] are odd under C. 

Given the ordering in eq. (|3.1|), the second quantized Hamiltonian becomes T~L[(f)} = | / dx 
[W, i?(0)\I/], where H(4>) is the single-particle Dirac Hamiltonian in the presence of <p(x). 
Denoting the (positive and negative) energy eigenvalues of H by u> n , and the corresponding 
eigenfunctions by ip ni we make the usual Fock decomposition, *&(x, 0) = >0 b n i[) n (x) + 
Hu n< Q d bl>n{x), so that 

H[$] = \ Un{blb n - b n b{) + \ Un(d n dl - d\d n ) . (3.2) 

OJn>0 W„<0 

We are interested in the expectation value of W[0 ] in states of definite occupation number. 
Consider the vacuum state |0), characterized by b n \Q) = d n \Q) = 0. Then 

E F {$] = (n\ni$]\n) = -lJ2\ UJ n\- ( 3 - 3 ) 

n 

Of course this sum is divergent. Its regularization and renormalization are discussed below. 

It is easy to confirm that \ is the lowest energy state. We see that the vacuum energy 
in a charge conjugation invariant theory is given by the (regularized) sum over the absolute 
values of all single-particle states. This result can be verified [H| using functional integral 
techniques, where it follows from the evaluation of the path integral for large Euclidean 
times. The fermion number of \Q) depends on 0, and for many configurations of interest it 
is not zero. 

For completeness, we write out the one-loop fermion vacuum contribution to the effective 
potential, 




Of course, our principal motivation in this work is to go beyond the effective potential 
approximation for the fermion one-loop contribution to the energy. Nevertheless, eq. ( |3.4|) 
will be useful to verify our numerical results for the fermion one-loop vacuum contribution 
in the case of a slowly varying boson fields. Also, like the boson effective potential, Vjg 
vanishes on the chiral circle, as we anticipated in the discussion of the a — > limit in the 
previous section. 
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B. Solutions to the Dirac Equation in a Chiral Background 

The Dirac Hamiltonian is 

H[$] = zo-i— + Ga 2 (f) 1 {x) + Ga 3 <p 2 (x)- (3.5) 

Although the underlying theory is charge conjugation invariant, the Dirac Hamiltonian in 
the presence of a fixed 4>(x) is not. Thus it is necessary to consider both positive and negative 
energy eigenvalues u of the time-independent Dirac equation 

Hip = uip. (3.6) 
The associated second-order equations for the Dirac spinor, ip = (f), are 

- /" - G</>'J + G^uj + G<j> 2 )-\f + G<j> x f) = (uj 2 - G 2 <t>l - G 2 <f> 2 2 )f 
-g" + Gfag - G^uj - Gfo)-\g ! - Gfag) = (to 2 - G 2 ^ - G 2 <$)g . (3.7) 

In one dimension, there are two channels for each energy. The S'-matrix is 2 dimen- 
sional and, in general, not diagonal. We simplify this situation using parity invariance. By 
demanding that 0i and <p 2 are respectively even and odd under coordinate reflection, i.e., 
<j>i(x) = <f>x(—x) and 4> 2 {x) = —4> 2 (—x), we ensure that the Dirac Hamiltonian, eq. (|3.5f) , is 
invariant under parity, 

[P, H] = 0, where P = 7o n = a 2 U (3.8) 

and n is the coordinate reflection operator that transforms x to —x. Thus we can choose a 
basis of parity eigenstates, 

PiP ± (x) = <J 2 iP±(-x) = ±*l> ± (x) . (3.9) 

Using parity we can replace the scattering problem on the line x G [— oo, oo] by two scat- 
tering problems on the half-line x G [0, oo] corresponding to even and odd parity. Eq. ( |3.9| ) 
gives boundary conditions at x = on the parity eigenstates: 

V>+(0) oc Q and V-(0) oc ^.J . (3.10) 

The solution to the Dirac equation on the half-line with either boundary condition is unique 
up to an overall normalization. For x — >• oo, this unique solution can be written as a 
superposition of incoming (oc e~ tkx ) and outgoing (oc e tkx ) waves. The coefficient of the 
outgoing wave relative to the incoming wave defines the phase shift. 

To implement this program, we introduce eigenstates of the free Dirac Hamiltonian with 
energy u, 



UJ 
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where k = +^/uj 2 — m 2 . Next, we construct the eigenstates of the full interacting Dirac 
Hamiltonian with energy u that are asymptotic to ip° ±k as x — > oo, 



/(*) 

,uo + G(j) 2 (x) 



V+k{x) - ( i {f{x) + GMx)f{x)) ) (3-12) 



and 



/*(*) 

= I {fix) + G(pi(x)f*(x)) I (:U:1) 

M + G<p 2 {x) 

where f(x) is the solution to the real second-order equation, eq. (|3.7|) , for the upper compo- 
nent, subject to the boundary condition that f(x) — > e lkx as x — > oo. It is easy to verify that 
in the same limit ip±k{x) — > <p± k (x) since the boson fields approach their vacuum values. 

Since the even and odd parity channels decouple, the S-matrix is diagonal. Its diagonal 
elements S± = e 2tS± ^ can be defined through the even and odd parity eigenstates of H, 

4>+(x) = f- k {x) H S + (oj)ip +k (x) 

UJ 

ip-(x) = tp-k(x) S-(u)ip +k (x) . (3.14) 

UJ 

If we set the interaction to zero (0i = v , 02 = 0), then if) + (■?/>_) reduces to the even (odd) 
parity solution to the free Dirac equation with S± — 1, which explains the extra factor of 
^ in eq. ([TTJ). 

To determine S± we use the fact that the eigenstates of eq. (|3.14j) obey eq. ( |3.10| ). For 
the positive parity channel this yields 

q(] _ (m+^) [(u-GMWM-rm 

+( J " .[(.-GMomw-m] (3 5) 

and similarly for the negative parity channel 

() " W [(a, + Gfc(0))/(0) + /'(0)] ■ (3 - 16) 

To compute the phase shifts efficiently and avoid 2ir ambiguities, it is convenient to 
factor the free solution out of f(x) by writing f(x) = e ikx e ^( x > w ) . Then the phase shifts are 
given by 

o±(v) = -Re/3(0,w) - arg H ■ — (3.17) 

L =pcj + m + ik J 

where the complex function /3(x,u) solves the differential equation 

- i/3"(x, u) + 2k/3'(x, uj) + (3' 2 {x, uj)-m 2 + G 2 <p\(x) + G 2 <g{x) - G<j>' x {x) 

+ Tnl) \ [G<f>i(^)+i(k + i3'(x,u J ))] = (3.18) 
u + G<p 2 {x) 

subject to the boundary conditions /3(oo,u) = P'(oo,u) = 0. It is this equation that we 
solve numerically for a given background 4>(x) to determine the phase shifts. 



8 



Heavy Fermion Stabilization of Solitons in 1+1 Dimensions 



C. Fermion One-Loop Contribution to the Effective Energy 

Having explained how to obtain the phase shifts for positive and negative energies and 
parities, we next show how to compute the fermion quantum corrections to the soliton energy 
following the method of 0-^. 

We are interested in the "effective energy" of a classical field configuration, which is the 
effective action per unit time for a time independent field, and is the energy of the lowest 
state with a given expectation value (<p(x)). For a self-coupled quantum field, as considered 



in Ref. [14], one introduces a source that makes the configuration <p(x) a stationary point of 
the classical action. A Legendre transformation then converts the energy from a functional 
of the source to a functional of the classical background field |[15| . For the case at hand, 
however, we have suppressed the quantum effects of the background field through the large- 
Np expansion. As a result, we can simply fix the classical backgound, at a cost in energy 
E c \ that is given by the classical Hamiltonian, and then compute the quantum corrections 
due to fermion fluctutatons in this background. 

The fermion effective energy in the presence of the classical background AE F [(p] is given 
by the sum over fermion loop diagrams with all sets of insertions of the background field. 
These diagrams can be evaluated by standard field theory techniques. Since some of them 
will diverge, we must introduce a regulator and counterterms that depend on the regulator. 
We will use dimensional regularization; other choices, such as a momentum cutoff, should 
be equivalent (in general, we expect two different regulators to lead to the same physical 
predictions as long as they preserve the same symmetries). The counterterms are then fixed 
by renormalization conditions on the parameters of the theory. Rather than sum all the 
diagrams, our approach will be to write an equivalent expression for the sum in terms of phase 
shifts, their Born approximations, and a small number of renormalized diagrams (which in 
this case will happen to be identically zero). This transformation relies on the correspondence 
between the Born expansion for the phase shifts and the diagramatic expansion of the 
effective energy. By comparing the two expressions for the energy as analytic functions of 
the regulator, we will verify that in making this correspondence, we have not introduced 
any finite errors or changed the definition of the theory. 

Formally, AE F [(f)] is also given by the shift in the zero-point energies of the fermion 
modes, — \ XXM ~~ l^ol), where {^o} are the eigenenergies in the free case. We rewrite this 
formal expression as a sum over bound states and an integral over k weighted by p{k) —p (k), 
the difference between the density of states in the free and the interacting cases, giving 

AE F [$] = - (j J2 K'l + \ J" dk ^+^(p(k) - p (k)) - f ) + A£ ct (3.19) 

where ui G [—m,m] denote the discrete bound states. The extra m/2 takes into account 
the contribution from the "half-bound" states at u — ±m in the free case |l|-||. AE ct 
is the contribution from the counterterm Lagrangian, which is evaluated at tree level. A 
counterterm proportional to (f)-(j) is implicitly present in eq. ( |2.1|) , which must be renormalized 
due to the coupling to the fermions. 

Both the counterterms and the integral over the continuum are formally infinite. Since 
our model is renormalizable, the sum is finite. In order to avoid ambiguities associated 
with the separate divergent terms in eq. ( ]3.19|) we have extended the method of dimensional 
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regularization to this representation of the effective energy: Eq. ( |3.f9| ) and the subsequent 
development should be understood as defined in n dimensions, where all integrals are con- 
vergent and all counterterms are finite. All the manipulations that we perform in the rest 
of this section are unambiguous in n dimensions. The final result is then analytically con- 
tinued back to 1+1 dimensions, where it remains finite. Since this extension of dimensional 
regularization is new and may have wider application, we present it in a self-contained way 
in Appendix [A]. 

The difference between the free and interacting density of states is given in terms of the 
scattering phase shifts by 

p (k) - p (k) = ^5 F (k) (3.20) 

where 

5 F (k) = 5 + (cu k ) + 6+(-u k ) + <L(w fe ) + 6-(-u h ) (3.21) 



and uj k = +\fk 2 + m 2 . 5p(k) sums the contributions from positive and negative energies 
and both parities. As in [0-0], we use Levinson's theorem 

+ --) 



S + (m) + 8 + (—m) = n{n 



<L(m) + &_(-m) = 7r(n" - ~) (3.22) 

to relate the phase shifts at k = in each channel to the number of bound states in that 
channel. The factors of ~ are peculiar to one dimension and reflect the existence of "half- 
bound" states at threshold in the absence of interactions. For a complete discussion, see 
Refs. and fl6|] . We express Levinson's theorem as 

, \ m f°° „ d „ m , 

in order to rewrite 

AE F [$] = - K-| + ±J~ dkVtfT^^6?{k) - y) + AE * ( 3 - 24 ) 



as 



AE F [$] = - (i (N - m) + jf °° ^ [VW+m 2 ' - m) ^{k)) +AE ct . (3.25) 
This substitution, which is valid in arbitrary dimensions, will allow us to formulate a Born 



expansion for eq. (|3.25|) without infrared divergences. 



In terms of the shifted field h = <pi — v, the <fi ■ <fi counterterm contribution to eq. Q3.25] ) is 



AE ct = C J( h2 + 2h v + <t>l) dx (3.26) 
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consistent with the chiral symmetry of the fermion interaction. We fix the coefficient C 
by the renormalization condition that the counterterm exactly cancels the tadpole graph 
(the term linear in h), and we perform no further finite renormalizations. This choice 
fixes the counterterm contribution to be equal to the tadpole graph plus those parts of 
the graphs with two external lines that are related to the tadpole graph by eq. (|3.26| ). 



This contribution is entirely local, i.e., independent of the external momenta. In order to 
reexpress this contribution in terms of phase shifts, we use the Born series. We expand the 
phase shifts in powers of the external fields h and 02- The contribution from the tadpole 
graph corresponds to the contribution from the first Born approximation to the phase shift. 
Expanding eq. ( |3.17| ) and eq. ( |3.18| ) to lowest order in h and 2 we find that the first Born 
approximation is given by 

6$\k) = + 5^\-iu k ) + S^(u k ) + (-<*>*) = — jf h(x) dx. (3.27) 

Once we have fixed the term linear in h, the entire counterterm contribution to the phase 
shift is fixed by eq. ( p.26|) to be 



S F (k) = / dx ( v 2 - 0(x) 2 ) . (3.28) 



Subtracting 5p{k) from 5p(k) in eq. ( |3.25|) then implements the full contribution of the 
counterterm, giving 

AE F [$] = -lJ2 (N " m ) ~ I" ^ (v^ 2 + m 2 - m ) ^ (<k(*0 - M*0) • (3-29) 



By expanding eq. (|3.17|) and eq. (|3.18| ) for large k we see that <5f(&) gives the leading 1/k 



behavior of 5p(k) for k large, and thus the resulting integral over the continuum is finite. 
We note that on the chiral circle 2 = v 2 , the counterterm contribution ( |3.28 ) to the energy 



vanishes, implying that the one-loop quantum contribution to the energy is finite. This is 
a consequence of the chiral invariance of the boson-fermion interaction, which forces the 
counterterms involving the "would-be" Goldstone boson to contain at least two derivative 
operators. On dimensional grounds these counterterms must be finite in 1+1 dimensions. 
In the present context, it implies that 5f(&) decreases more rapidly than 1/k as k — > oo for 
scalar configurations on the chiral circle, and our numerical computation indicates that in 
this case 5f(&) goes like 1/k 3 for k large. 

For numerical computations, it is convenient to integrate eq. (|3.29|) by parts, yielding 

^ r 00 dk k 

A£ F [0] = 4E(N-™) + / o ^ VWT - ? (S F (k) - S F (k)) . (3.30) 

This is the finite expression for the fermion one-loop contribution to the effective energy 
which we evaluate numerically. The analysis of the first part of this section gives us the 
phase shifts, and we find the bound states by ordinary shooting methods, with Levinson's 
theorem telling us how many to look for. 
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D. The Fermion Number of a Configuration 

We are interested in fermionic solitons and therefore must keep track of the fermion 
number. If, as we adiabatically turn on the background configuration starting from the trivial 
vacuum, a fermion level crosses zero from above, the vacuum will have this state filled, giving 
the configuration a positive fermion number. (If the crossing is in the other direction, there 
is a negative fermion number.) The fermion energy levels are iV^-fold degenerate, so the 
configuration carries each of the Np fermion numbers. We know from general considerations 
that if <f>(x) circles = (0, 0) as goes from (v, 0) at x — — oo to (v, 0) at x — oo, then the 
vacuum state will carry nonzero fermion number provided that w, the scale over which 
varies, is much larger than the fermion Compton wavelength 1/m HjlTj • 

In Ref. 0] we show how to compute this fermion number directly using Levinson's the- 
orem, which relates the phase shifts at u — ±m to the total number of bound states. We 
use that the difference between the phase shift at the positive energy threshold, u = m, and 
the actual number of positive (negative) energy bound states n> (n<) counts the fermion 
number Q of the vacuum: 

Q\4>\ =N F (^(8 + (m) + 54m)) +§-n>) 

= -N F (-(6+{-m) + <L(-m)) + \ -n<) . (3.31) 

As noted above, the extra | is peculiar to one dimension, and is essential to get correct 
(integer) values for Q. "Half-bound" states at the respective thresholds count as | in ri> or 
n < . The configurations that we consider circle = (0,0) at most once, so Q is either or 
Np. 

We are interested in states with fermion number Np. If Q = Np, then the state we want 
is the vacuum. If Q — 0, then we build the lowest energy state with fermion number Np 
by filling the lowest positive energy level uj\ of eq. ( |3.5| ) with Np fermions. Thus the total 
fermionic contribution to the energy is 

E F = AE F + uj x (N F - Q) . (3.32) 

We note that E F varies smoothly with uj\ even as iO\ crosses zero. The vacuum and 
valence contributions each have discontinuities in slope when u>i = 0, but the contribution 
from this level is — \\uj\ \ from the vacuum energy in eq. (|3.30| ) plus ui9(uji) as a filled valence 
level. The sum is |u;i for all u)\, which is smooth. 

The spectrum also contains states with total charge between zero and Np. These states 
are constructed by an appropriate filling of the tightly bound level if u)\ > or the emptying 
of that level if uji < 0. We do not consider such states because their binding will be smaller 
than that of states with charge A 7 ^. We note, however, that if the effects we find in the large 
Np limit persist for moderate values of Np, then in that case a stable soliton with fermion 
number 1 would represent an elementary fermion. 

It will also be helpful to consider the fermion charge density. The vacuum contribution 



is given by summing over the eigenstates obtained from eq. ( 3.14 ). 



(n\^(x)^(x)\n) = -N F I sgn(cj)±(V+(a;)V^) + ^-(^)V-(^)) du (3.33) 

J — oo 
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where the integral over u includes both continuum and bound states. We have restored the 
suppressed u label on the wavefunctions and normalized them by 

il>Z(x)ty*[(x)dx = 5(w-u/) (3.34) 

where the right hand side is interpreted as a Dirac delta function for continuum states and a 
Kronecker delta for bound states. The norm of the scattering states is related to the phase 
shifts by 

Ul{x) ty±(z) -l)dx = I^M (for | w | > m \ (3.35) 

71 duj 

as derived in [TJJ. The spatial integral over the bound state contribution to fl3.33| ) yields the 
difference between the numbers of positive and negative energy bound states, n > —n < . Thus 
by integrating eq. ( p. 33 ) over x and using Levinson's theorem and the fact that the total 



number of states is unchanged, we obtain eq. ( |3.31| ) 



If Q = we have to explicitly add the contribution of the level ipi(x) with the lowest 
positive eigenvalue. Thus the total charge density is given by 

j (x) = (Np - Q) ifil(x)ik(x) + (Sl\tf{x)*(x)\n) . (3.36) 

4. THE TOTAL ENERGY 

Now that the theory has been prescribed by Lb + £f, it is useful to introduce dimension- 
less variables to simplify our analysis. We measure all energies in units of the perturbative 
fermion mass, m = Gv, and all distances in terms of £ = mx. We are interested in calcu- 
lating the energy of a fixed background configuration <p(x), which we describe in terms of 
a small set of dimensionless variational parameters {0}- Next we rescale by v (both of 
which are dimensionless in one dimension) so that <fi — > (1,0) as |£| — ► ±oo, and we define 
dimensionless couplings 

61 = aIld X = ^2 ' ( 4 - X ) 

With all these rescalings the boson classical energy can be written as 



V 

m 



= v 2 £ d (a, A, {Q}) (4.2) 

where prime denotes a derivative with respect to £. 

The fermion one-loop contribution to the energy is given by a regulated sum over the 
absolute values of the eigenfrequencies, {uj n }, of H[(f>}. Once we measure distances in units 
of £ and choose dimensionless variational parameters {Q}, the perturbative fermion mass 
m factors out of the Dirac Hamiltonian, so its eigenvalues {uj n } scale with m. In all, the 
fermion one-loop vacuum energy can be written as 

E F $] = N F m£ F {{C l }). (4.3) 
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Thus we compute 



E tot \<t>] v 



2 



— £ d (a, A, ({&})• (4.4) 

JSp m JSp 

This computation is exact in the limit Np — > oo and v 2 — > oo with the ratio v 2 /Np held 
fixed, since the contributions we have neglected all arise from loops with internal bosons, 
which are suppressed by v 2 relative to the classical energy and by Np relative to the one-loop 
fermion energy. 

To compare the energy of our background configuration to the energy of the state with 
the same charge built on top of the translationally invariant vacuum, we study the effective 
energy per fermion minus Np times the perturbative fermion mass Gv, which we denote as 
B[4>\. In terms of rescaled variables, 

b(^, a, A, {0}) = ^ci(«, A, {0}) + £ F ({Q}) - 1 • (4.5) 

A field configuration with B < is energetically favored over the state with the same charge 
Np built on top of the translationally invariant vacuum. 



5. NUMERICAL EXPLORATIONS 

In this section we present the results of our numerical studies. We will show that the 
model has a stable fermionic soliton for a wide range of the model parameters A, a, and 
v 2 /Np. For fixed model parameters, we search over the variational parameters looking for a 
bound solitonic fermion. The classical contribution to the soliton energy, £ c \ in eq. ( |4.5|) , is 
positive and is scaled by v 2 /Np. The fermion one-loop contribution, S F — 1, is generically 
negative. Thus, the existence of a stable fermionic soliton can be discussed in terms of a 
maximum value of v 2 /Np for a given choice of A and a. 

A. Variational Ansatz 

— * 

We choose an ansatz <fii(£) characterized by the parameters R and w, 

fafaRtW) = (1 - J R + J Rcose / (^,w), J Rsine / (^,w)) (5.1) 

with 

ej(f,w) = 7r(l + tanh(f/u;)). (5.2) 

For fixed R and w, (pj describes a circle centered at 1 — R with radius R as £ varies, w gives 
the characteristic size of the configuration. Varying R allows us to interpolate between the 
trivial vacuum, = (1,0), and a configuration that circles the origin once while staying on 
the chiral circle, \<f>\ = 1. Configurations with R > 1 go beyond the chiral circle. When the 
symmetry breaking parameter a is large, <fi = (—1,0) is energetically disfavored, leading to 
a preference for R to go to zero to avoid that point. 

We have also considered other ansatze, but we found that they yielded variational minima 
that had comparable or higher energy. Thus we will restrict our attention to this choice. 
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B. Stability Studies 



To search for the soliton, we must compute the classical energy £ c i(oc, X, R, w) and the 
fermionic contribution £ F (R, w). Our scaling has simplified the numerical ansatz by putting 
all the dependence on the model parameters A, a, and v 2 /Np into the classical energy, which 
is easy to compute. The fermionic contribution, which is harder to compute, depends only 
on the variational parameters R and w. 

The classical energy is computed from eq. Q4.2|) . A typical energy surface for ansatz <pi is 
shown in Fig. [I] as a function of R and w for a generic choice of model parameters, A and a. 





FIG. 1. The classical energy, £ c \, as a function of the variational parameters R and w in eq. (4.2 ) 
The parameters of the hosonic Lagrangian are X = 1.0 and a = 0.5. 



Simple scaling arguments allow us to understand the principal features of Fig. |l[ For fixed 
w, £ c i(<5, A, R, w) vanishes as R — ► because the ansatz approaches the classical vacuum for 
all £. For fixed R ^ 0, £ c \(a, A, R, w) diverges like 1/w for small w due to the |0'| 2 terms in 
eq. ( fOl) . For large w, £ c i(oc, A, R, w) grows linearly with w for fixed R ^ 1 because V(<f>) is 
nonzero in a region of size w. 

Next we turn to the fermion contribution to the energy. There are two distinct contri- 
butions discussed in the previous section: first, the energy of any filled "valence" level, and 
second the vacuum energy. An important property of our ansatz is its propensity to bind a 
fermion. This can be seen in Fig. where we plot the lowest fermion eigenenergy uj\ as a 
function of the variational parameters R and w for the ansatz <fri. Since we are using scaled 
variables the fermion continuum begins at uj = 1. Notice that for R near 1, which puts <pi 
close to the chiral circle, the lowest fermion eigenenergy decreases quickly with increasing 
w. Even a modest value of w leads to a negative eigenenergy. 

The strongly bound fermion state drives the formation of a fermionic soliton. The situ- 
ation is reminiscent of the appearance of a fermion zero mode in the presence a soliton in a 
theory of a real scalar field |T8 |. It is expected from considerations of level crossings in the 
presence of an adiabatically changing field: If we fix R > ^ so that the orbit in the <p2) 
plane circles the origin, then as w increases from to oo one fermion level must cross from 
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FIG. 2. The lowest quark eigenenergy, uj\, for ansatz as a function of R and w. Note that 
for large R and w, uj\ is negative. A solid curve marks the contour oj\ = 0. 



the positive to the negative continuum |6|,|n]]. Our numerical studies indicate that iO\ is 
never negative unless the orbit of (f) circles the origin. 




0.2 

FIG. 3. The vacuum contribution to the one-loop fermion energy as a function of R and w. 
The sharp edge in this contribution occurs when the lowest fermion eigenvalue crosses zero. See 
Fig. § 

In Appendix [FJ, we explore the generality of level crossings for field configurations on the 
chiral circle. We show that the WKB approximation yields a single strongly bound level in 
the positive parity channel, which crosses from the positive to negative continuum as the 
width is increased. Note that merely examining the spectrum in the limit of large width 
gives no evidence of the level crossing. It is essential to follow the levels as a function of 
width or use the Levinson's theorem method, eq. (|3.31| ). 

Next we calculate the contribution from the fermionic vacuum - the Casimir energy - 
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given by eq. ( ^.30[ ). Note that this sum over the eigenenergies of the Dirac Hamiltonian is 
of the same order in v 2 /Np as the energy of the lowest eigenmode and must be included in 
a consistent calculation of the one-loop fermion contribution to the effective energy. This 
piece of the energy is also a function only of the ansatz parameters R and w. It is shown in 
Fig. ||. It displays a discontinuity in slope at the point where the lowest eigenvalue passes 
through zero. Fig. shows that the sum of the Casimir energy and the contribution from 




FIG. 4. The fermion one-loop contribution, £ F (R, w), to the energy of a state of fermion number 
Np for the ansatz 4>j. 



uj\ is smooth, in agreement with the discussion in the previous section. 




FIG. 5. B as a function of the ansatz parameters for a = 0.5, A = 1.0, and v/y/Np = 0.375. 
A solid curve marks the contour B = 0, and a star indicates the minimum at w = 2.808 and 
R = 0.586. 

Finally, we combine the classical energy and the one-loop fermion energy to form 
B(v 2 /Np, <i, A, {0}) and search for a stable fermionic soliton. As an example of our re- 
sults, Fig. H shows B as a function of R and w for A = 1, a = 0.50, and v/y/Np = 0.375. 
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FIG. 6. (pi, (f>2, and the fermion number density jo at the variational minimum for 5 = 0.5, 
A = 1.0, and v/^/I^ = 0.375, which is at R = 0.586, w = 2.808. The left panel shows 
and 4>2(0 at this point, and the right panel shows jo(0> which gets contributions from both the 
fermion vacuum and the filled valence level. 



B < is the signal of a stable fermionic soliton. The domain of negative B is large and 
the variational approximation to the fermionic soliton lies at the minimum: R = 0.586, 
w = 2.808. In this case the lowest fermion eigenenergy, uj\ = 0.0985, is slightly positive 
so this level must be filled. The variational approximation to the binding energy is 0.253. 
In Fig. H we show <fii(£), faiC), and jo(x). We see that the charge density is well localized 
around the center of the soliton. 

We are now in a position to study the stability of the fermionic soliton as a function 
of the model parameters A, a, and v 2 /Np. For each choice of parameters we minimize B 
over the ansatz parameters, R and w. The output is B(X,a,v 2 /Np) and the parameters 
R(X, a,v 2 /N F ) and w(\, a,v 2 /N F ) that minimize B. In Fig. ^ we plot B versus v/s/Np for 
various choices of A and a and see that binding occurs over a wide range of model parameters. 
Binding is a quantum phenomenon - it is maximal when v/y/Np, which multiplies the 
classical contribution, vanishes, and decreases with increasing v/y/Np. 
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FIG. 7. B as a function ofvj\J Np for various values A with a = 0.25 (left panel) and for various 
values a with A = 1.0 (right panel). 
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FIG. 8. The width ('left panel) and the radius (right panel) of the configurations that minimize 
the total energy as a function of the explicit symmetry breaking a. Several values of the Higgs 
coupling constant A are considered and v/y/Np = 0.375. 



The dependence on the explicit symmetry breaking parameter a is particularly important 
because of the infrared problems we expect to arise as a — > 0. Fig. [7| confirms our expectation 
that B — > — 1 (strong binding) as a — > 0. However this is a suspect limit. A look at the 
width of the soliton displays the problem. This is shown in Fig. || where w and R are plotted 
versus a for various choices of A and a typical value of v/y/Np. Two qualitatively different 
regimes are separated at a ~ 0.3. Above a the size of the soliton is almost independent 
of a. Below 5lq the width of the soliton diverges like l/y/di and the ansatz moves towards 
the chiral circle, R — l. However, this is where we expect the model to be invalidated by 
infrared divergences generated by the would-be Goldstone mode. Therefore we only trust 
our results in the domain a > a - We have made plots analogous to Fig. [| for various 
choices of A and v/ y/N F . We find that <5 depends weakly on A but strongly on v/ y/N F . In 
Fig. [5] we plot a (v/y/N F ) defined as the point where the break occurs in figures analogous 
to Fig. p. In the absence of a complete analysis of the infrared problem as a —>■ 0, we only 
trust our results when a > 5>o{v /y/Np). 

For fixed a, when v/ y/Np is too small we run into the infrared instabilities we described 
above. When v/y/Np is very large the soliton binding becomes small. The interesting 
domain is between these two extremes where a stable soliton can be reliably described. To 
indicate this domain we have added a contour corresponding to 5% binding in Fig. || In 
between we have found a strongly bound fermionic soliton. 



6. DISCUSSION AND CONCLUSIONS 



We have shown that quantum effects can stabilize a soliton in a theory with no soliton at 
the classical level. In order to convincingly demonstrate this phenomenon, it is necessary to 
consider a renormalizable model. We study a (l+l)-dimensional model with a chirally in- 
variant Yukawa interaction between fermions and scalars. The scalar self-interaction is given 
by a Higgs potential with an explicit symmetry breaking term. (Without the explicit sym- 
metry breaking, the fluctuations of the Goldstone modes would restore the spontaneously 
broken symmetry.) We define the renormalized energy functional using counterterms de- 
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FIG. 9. The regions of soliton stability in the plane of v/y/Np and a. In the shaded area on 
the left, a growing width indicates potential infrared instabilities. In the shaded area on the right, 
the soliton is bound by less than 5 percent. In between, we have a stable, tightly bound soliton. 

termined in the perturbative sector of the theory. Our phase shift formalism allows us to 
unambiguously include the exact fermion one-loop quantum corrections in a manner that 
is computationally tractable. The omission of the scalar loop contribution is justified by 
assuming a large number of fermion flavors. 

Holding the total fermion number fixed, we performed a variational calculation to search 
for configurations with lower energy than the same number of free fermions. Once we have 
found a bound configuration we are assured that a stable soliton exists since its energy can 
only be smaller than the variational minimum. As expected, when the typical extension of 
the background field is large compared to the Compton wavelength of the noninteracting 
fermion, the fermion number is carried intrinsically by the background field, instead of 
through explicitly filled levels. 

We have found a wide range of model parameters for which we believe that the soliton is 
not destroyed by the infrared singularities of the "would-be" Goldstone mode, and where the 
binding of the variational minimum is sizable. In this case, the cost in the classical energy 
is more than compensated by the gain from the fermion loop. For a fixed set of the dimen- 
sionless model parameters, the gain in energy is proportional to the dimensionful Yukawa 
coupling constant G, and the phenomenon of soliton formation becomes more pronounced 
as the perturbative mass of the fermion increases. Thus, in 1+1 dimensions, heavy fermions 
can indeed stabilize solitons. 
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APPENDIX A: DIMENSIONAL REGULARIZATION 
OF THE CASIMIR ENERGY 

In the body of the paper we introduced and manipulated formally divergent expressions 
such as eq. (|3.19| ) without prescribing a regularization scheme to render them finite. For 
Feynman diagrams, regularization and renormalization are so familiar that such a casual 
approach would not warrant further comment. In this Appendix we show how to extend the 
method of dimensional regularization to Casimir sums, thereby rendering our earlier ma- 
nipulations well defined. This analysis resolves ambiguities in the definition of the effective 
energy that have been noted in recent works [0. 

We begin by defining the Casimir energy in n (space) dimensions in terms of a sum over 
bound states and an integral over a continuum density of states. These are properties of 
a radial Schrodinger-like equation, which remain well defined when n is noninteger. We 
choose n such that the expressions are finite. We then isolate and compute the term AE^ 
that diverges when n is continued to n = 1, the case of interest. For < n < 1 we 
show explicitly that AE^ is identical to the contribution of the lowest order Feynman 
diagram in the loop expansion of the effective energy. This is our fundamental result. Still 
keeping < n < 1 we subtract AE^ l > from the Casimir sum and add back in the Feynman 
diagram. The subtracted Casimir sum is finite for n = 1 while the now (conventionally) 
dimensionally regulated Feynman diagram is renormalized by the counterterm contribution, 
AE ct . This calculation demonstrates that we have precisely implemented the standard 
dimensional regularization and renormalization process of perturbative field theory. 

For simplicity we consider the self-interactions of a single real boson in one dimension; 
the generalization to fermions is discussed in Ref. HI]. We take the bosonic Lagrangian 

c = ±dt l( f>dy-u( ( p) + cu"((j)) (A.i) 

where prime denotes differentiation with respect to 0. We have indicated the counterterm 
CU"((f)) explicitly including its cutoff-dependent coefficient C, and we are considering an 
arbitrary potential U(<p). This counterterm renormalizes the boson's mass, and is the only 
counterterm needed to render the theory finite in one spatial dimension. We take a back- 
ground (po(x) that is either a solution to the classical equations of motion or held in place 
by an external source, so that it is a stationary point of the classical action. We wish to 
compute the one-loop contribution to the energy arising from the bosonic fluctuations about 
4>o(x). To simplify the analysis, we assume it is either an odd or an even function of x, and 
we take U(<p) to be an even function of <ft that gives a perturbative mass m to <fi. The 
potential for small oscillations around <fio(x) is then given by 

V(x) = U"(Mx)) ■ (A.2) 

Our restrictions on U((f)) and (f>o(x) ensure that V(x) is even in x. As usual, the relationship 
between the scattering phase shifts and the density of states is 

Sp(k) = p{k) - Po (k) = — Tr lnS = . (A.3) 

For a single real field 0, we only need to consider positive energies. We use parity to 
decompose the 5-matrix into symmetric and antisymmetric channels, giving 

5 B (k) = 6 + (k) + S-(k). (A.4) 
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As in Section III we define the Casimir energy as the sum over bound states {ojj} plus an 
integral over the continuum weighted by u = \/k 2 + m 2 , 



POO 

AE[$) = \ - m) + \ \ dk{u - m)Sp(k) + AE ct 

Jo 



(A.5) 



where in analogy to eq. ( |3.19| ) we have used Levinson's theorem to regulate potential infrared 
divergences (as k — > 0). 

We now generalize eqs. (|A.3|) and ( A.5|) to n dimensions. The scattering problem gener- 
alizes to a central potential in n dimensions. The S'-matrix is diagonal in the basis of the 
irreducible tensor representations of SO(n). These are the traceless symmetric tensors of 
rank £, where £ = 0, 1, 2, ... . The formula for the density of states, eq. ( |A.3j ), becomes a 
sum over £, 

-i j i i 00 

SMk) m = j2 NXt(k) (A . 6) 

e=o 

where 8 n ,i{k) is the phase shift in the £ th partial wave and is the degeneracy of the SO(n) 
representation labelled by £. For integer n and £, is given by the dimension of the space 
of symmetric tensors with £ indices that each run from 1 to n, with all traces (contractions) 
removed. Working out the combinatorics gives 

= ( n + £-l)\ _ (n + l-3)! 

£\{n-l)\ {£-2)\{n-l)\ ' 1 ' ' 

To prepare the way to continue to noninteger n we write N„ in terms of T-functions instead 
of factorials, 

n = rKni+i) < n+ "- 2) - (A ' 8) 

For n = 3, iV„ reduces to 2£ + 1 as expected. 

The phase shifts are obtained by solving the radial Schrodinger equation generalized to 
n dimensions, 

- ^—Itf + ^ + n ~ 2 \ + 2mV(rU = k 2 ip (A.9) 

which is related to Bessel's equation for V = 0. At the origin, the regular solution ip n ^ is 
proportional to r independent of n. 

Incoming and outgoing waves are generalizations of spherical Hankel functions, 

h ( nf(kr) = — -L-j (J^_i(*r) ± iYn + ^{kr)) (A.10) 
[Kr) 2 

and the phase shifts are defined in the usual way by writing the solution ip n ^ regular at the 
origin as 



~ h%\(kr) + e 2iS ^hW(kr) (A.ll) 
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for large r, where the potential vanishes. The leading behavior of S n> e{k) at large k is given 
by the first Born approximation, 



Jn + £-i(kr) 2 V (r)r dr . 



(A.12) 



The expression for the Casimir energy in n dimensions is 



oo 

AS n [0] = iE(E iV n< 

j e=o 



CO 



m 



— (u-m^NZ—SnM+AE^ (A.13) 
71 1=0 



where N* is given by eq. ( |A.S| ), the 0Jj, n ,i are the normalizable solutions to eq. ( |A.9| ) in 
each partial wave £, and 8 n ,t{k) is determined from eqs. ( |A.9|) and eq. (|A.11[) . AE ct ^ n is the 
counterterm contribution, to be fixed by a renormalization condition below. The rest of 
eq. (|A.13 ) is well defined for noninteger n, where the integration over k and the sum over £ 
converge. Holding the divergences in abeyance we can verify that eq. ( |A.13| ) reduces to the 
naive result, eq. ( |A.5| ), as n — > 1: N[ vanishes for all I except I = and £ — 1, where it is 1; 
for £ = 0, 5i y o(k) and u>j,io are obtained from the solutions to eq. ( |A.8| ) that have vanishing 
first derivative at r = 0, while for £ — 1, <5i,i(fc) and i i are obtained from the solutions to 
eq. ( |A.^ ) that vanish at r = 0. Thus £ = and £ = 1 correspond to the even and odd parity 
channels respectively. 

Our approach consists of subtracting the first Born approximation to the phase shift and 
replacing its contribution by the tadpole graph, which we then calculate in ordinary Feynman 
perturbation theory. Thus we must demonstrate explicitly that these two quantities are equal 
by computing both as analytic functions of n, away from integer n, where both diverge. Once 
the leading Born approximation has been subtracted, the integral over the phase shift is finite 
and we can take the limit n — > 1 with no further subtleties. 

The tadpole graph requires the external momentum to be equal to zero. Thus we should 
expect that both the leading Born approximation and the tadpole graph will depend only 
on the spatial average of the potential, 



(V) 



V{x) d n x 



27rf 

r(f) 



V{r)r n ^dr. 



(A.14) 



From eq. (|A.13|) , the contribution to the energy from the first Born approximation is 



d5^(k) 



dk 



(A.15) 



We obtained (u — m) rather than uj in the integrand of eq. (|A.15| ) because of our use of 
Levinson's theorem to put the Casimir energy in the form of eq. (|A.5|) . This manipulation 
was required to ensure the absence of infrared singularities in one spatial dimension, but we 
will see that it is also necessary for us to write a sensible expression in arbitrary dimensions.^ 



the boson were massless, we would have found an infrared divergence in eq. (|A.5| ) as n — > 1 from 
the 1/k divergence of the Born approximation at k = 0. This divergence reflects the well-known 
infrared divergences of massless theories in 1+1 dimensions [12]. 
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Using the Bessel function identity 



y (2g + 2l)r(2g + l) ()2 _ T(2q+l) ,z\* 



T(i+1) ™' T{ q +lf\2 

and setting q = ~ — 1, we explicitly sum over £ in eq. ( |A.15| ) and obtain 

AEW[$] = - ^ (n-2) [ X \ u - m)*"" 8 dk . (A.17) 

The k integral can be calculated in the vicinity of n = | and then analytically continued, 
yielding 

<._ mr -3^_^ r( I^ )r( ^). (A , 8) 



Hence we find 



(4tt)^ V 2 

which is exactly what we obtain using standard dimensional regularization for the tadpole 
diagram in n + 1 space-time dimensions. 

We choose renormalization conditions such that the tadpole graph is exactly cancelled 
by the counterterm contribution. Thus we can implement the contribution of AE ct ^ n by 
subtracting 5^\(k) from 5 n ^(k) in eq. ( |A.13| ), yielding a finite result 

OO /*00 J 7 < ~^~ > 7 

A^0] = §££*i("w-™)+ / ^("-™)£*%(M*0-tfj(*)) • 

(A.20) 

This result can then be smoothly continued to n = 1, giving 

A^[0] = I £fa - m) + j( -(w - ro) J- - *«(*)) (A.21) 

where 

= — / V(r)dr. (A.22) 

By continuing to fractional dimensions, we have regulated the theory, rendering it finite. 
We made certain that we held the physical renormalization conditions fixed while removing 
the regulator. This process defines the theory in terms of physical parameters that can be 
measured within the perturbative sector. Based only on these inputs, we can then calculate 
the energy of nontrivial field configurations. 
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APPENDIX B: STUDIES OF THE TIGHTLY BOUND FERMION LEVEL 

Our soliton's stability was driven largely by the strong binding of a single fermion level 
in a chiral background. In this appendix we explore this phenomenon further by studying 
the case of a scalar field constrained to the chiral circle, = (cos 6 (£), sin 6 (£)), with the 
chiral angle 0(£) varying monotonically between and 2tt as £ goes from — oo to +oo. We 
assume 0(£) = n — 6(— £), so 0(0) = tt and 6'(0) > 0. Using the WKB approximation, we 
argue that the vacuum will acquire unit fermion number when 0(£) varies slowly compared 
with m. 

It is convenient to adopt a basis for the Dirac matrices different from the body of the 
paper: 7 = o"i, 71 = io<i and hence 75 = 03. Denoting the upper and lower components of 
the spinor by a and b, respectively, the first order Dirac equations read 



-iB 



<i — e (cub + ib') 
e l ~ (ua — ia) 



(B.l) 



where as before, the energy u is measured in units of the fermion mass, m = Gv and a 
prime denotes a derivative with respect to £ = mx. In this basis, the second order equations 
become 



— a 
-b" 



(lu 2 -l)a + Q' (ia! -ua) 
(u 2 - 1)6-6' {ib' + ub) 



A bound state is a solution to eq. ( B.l ) which falls exponentially as £ 
parity: o(0) = ±6(0). 



(B.2) 
(B.3) 

00 and has definite 



Positive Parity 



Negative Parity 





1/0' 



FIG. 10. The bound state energies of the simple model as functions of the inverse of the constant 



slope 0' . In the positive parity channel the WKB solution (B.6) is also shown 



First, we suppose G is slowly varying and use the WKB approximation. We parameterize 
a(£) = e % f(€) and neglect f" compared with f' 2 . Substituting into eq. (|B.2|) and solving for 
/' we obtain 



0' 



2 



(B.4) 
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where we have chosen the root that gives a wavefunction that falls exponentially at large £. 
&(£) is given by the second of eqs. ( |B.1|) . The eigenvalue condition for positive (negative) 
parity at £ = reduces to 

-u - /'(0) = ±1 . (B.5) 

This equation has only one solution when 0'(O) > 0, 

w = -l + e'(0)/2 (B.6) 

which has positive parity. 1/0'(O) measures the width of the scalar field configuration. For 
large width, the WKB solution lies just above the negative energy continuum. As the width 
decreases, the WKB bound state energy increases, crosses zero when 0'(O) = 2 and enters 
the positive energy continuum when 0'(O) = 4. 

The condition that /" < f' 2 limits the validity of the WKB approximation to 8" < 1. 
For example, when uj = the condition reduces to 



0"(£) 



X 4-0' 2 (O 



< 1 . (B.7) 



This condition can be satisfied for all £ by making 0" small enough. Note 0"(O) = by 
symmetry thereby avoiding the apparent singularity at £ = where u = requires 0'(O) = 2. 

We have augmented the WKB analysis by solving a simple toy model where much of the 
calculation can be done analytically. In this model we let 0' = 7 =constant for — vr/7 < 
£ < 71-/7 and zero elsewhere. The positive and negative parity eigenstates are determined by 
a simple transcendental equation which can be solved numerically. The results are shown 



along with the WKB estimate, eq. ( |B.6|) , in Fig. [10]. As expected, the toy model has 
a positive parity bound state that descends rapidly from the positive energy threshold, 
through zero energy, toward the negative energy continuum. It is well approximated by the 
WKB estimate. The model has other positive and negative energy bound states, which are 
missed by the WKB approximation, but which remain in the vicinity of either the positive 
or negative energy continuum for all values of 7. 



APPENDIX C: REFLECTION COEFFICIENTS 

Ref. H claimed that the reflection coefficient for fermions scattering off a soliton back- 
ground should vanish. Since we calculate the phase shifts 5±(u), we can easily compute the 
reflection coefficient 

r(w) = \ [S+(u) - S-(u>)] (C.l) 

in terms of the S-matrix elements S±(u) = exp(2i5±(u)). We are thus equipped to reexamine 
this statement for our variational approximation to the soliton. In Fig. [11] we display the 
complex reflection coefficients as functions of the momentum of the scattering fermion for 
positive and negative energies. Although these coefficients approach zero very quickly as k 
increases, they do not vanish identically. In particular we find that these coefficients acquire 
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more structure as we get closer to the (unstable) symmetric formulation (5 = 0). The added 
structure is mainly caused by the increasing number of bound states in the wider potential. 
Because of Levinson's theorem these bound states cause the reflection coefficient to circle 
around the origin in the complex plane. Though it is still possible that the true soliton gives 
a reflectionless potential, we see no indication of this behavior. 



Re(r), (B > 
Im(r), co > 
Re(r), CO < 
Im(r), co< 



Re(r), 00 > 
Im(r), co > 
Re(r), CO < 
Im(r), co < 



0.50 

k 



0.50 

k 



FIG. 11. The reflection coefficients of the variational approximation to the soliton as functions 
of the momentum k = \/u> 2 — 1. We consider two values of a: 0.1 (left panel) and 0.6 (right panel). 
In both cases we have used A = 1.0 and v/y/Np = 0.375, as in Fig. ||. 

The statement in Ref. || is based on the correct requirement that the fermion current 
j\x) = (N F - Q)^i(x) f 7°7Vi(^) + (aOVV* (x)|ft) (C.2) 

must vanish for all x in the soliton background. However, decomposing \l/(x) in terms 
of creation and annihilation operators that multiply the solutions to the Dirac equation 
given in eq. ( |3.14|) , we find that j 1 (x) vanishes mode by mode as a simple consequence of 
\S + \ = \S-\ = 1 which is enforced by unitarity. Thus the requirement j l (x) = is trivially 
fulfilled and places no restriction on the soliton configuration. 

Though two of the most famous examples of topological solitons, the kink and the sine- 
Gordon soliton, both correspond to reflectionless potentials, we have no reason to assume 
that such a restriction applies in general. Indeed, the Jacobi model studied in Ref. |19| 
gives a family of solitons labelled by the elliptic parameter v. These solitons correspond to 
potentials that are not reflectionless except in the limit v — > and v — > 1 (where the model 
reduces to the sine-Gordon case). 

For certain interactions, the inverse scattering method applies and the equations of mo- 
tion (expressed in terms of scattering data) indeed yields reflectionless potentials pOfl . How- 
ever, the attempts to extend this method to chiral models in Ref. ||,10] are problematic: 
First, they cannot include the symmetry breaking terms needed to avoid the restoration of 
chiral symmetry by boson loops. Second, because these solitons approach different configu- 
rations as x —>■ — oo and x — > +oo (which are degenerate vacua in the absence of symmetry 
breaking), the methods developed in Ref. PJT7H show that the polarized Dirac sea carries 
fractional charge for each fermion flavor even though no level has crossed zero. When this 
contribution is properly included, the total fermion number of the configuration constructed 
in Ref. IHI adds to zero. 
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